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1 INTRODUCTION 



Both Galactic and extragalactic accretion-powered X-ray 
sources display significant aperiodic variability on a wide 
range of timescales. Similarities in behaviour between the 
black hole candidate X-ray binaries (BHCs) and the nu- 
clei of Seyfert galaxies (that presumably also harbour cen- 
tral black holes) have been used to suggest that the phys- 
ical processes occuring in accretion flows around the black 
holes in these objects are essentially the same, but scaled 
in some appropriate manner. There are indeed suggestions 
that such scaling relationships might be used as an estimate 
for black hole masses (e.g. Vaughan, Fabian & Nandra, 2002; 
Markowitz et al. 2003). Although X-ray binaries containing 
neutron stars as the accreting compact objects also exhibit 
interesting aperiodic, and quasi-periodic, behaviour (see for 
example the review by van der Klis 1989), we restrict our 
attention in this paper to black hole accretion discs, partly 
in the hope that understanding the physics of such accretion 
flows might be somewhat less complicated, but also because 



ABSTRACT 

Observations of accreting systems often show significant variability (10-20 percent 
of accretion luminosity) on timescales much longer than expected for the disc regions 
releasing most of the luminosity. We propose an explicit physical model for disc vari- 
ability, consistent with Lyubarskii's (1997) general scheme for solving this problem. 
We suggest that local dynamo processes can affect the evolution of an accretion disc 
by driving angular momentum loss in the form of an outflow (a wind or jet). We 
model the dynamo as a small-scale stochastic phenomenon, operating on roughly the 
local dynamical timescale. We argue that large-scale outflow can only occur when 
the small-scale random processes in neighbouring disc annuli give rise by chance to a 
coherent large-scale magnetic field. This occurs on much longer timescales, and causes 
a bright large-amplitude flare as a wide range of disc radii evolve in a coherent fash- 
ion. Most of the time, dynamo action instead produces small-amplitude flickering. We 
reproduce power spectra similar to those observed, including a 1// power spectrum 
below a break frequency given by the magnetic alignment timescale at the inner disc 
edge. However the relation between the black hole mass and the value of the break 
frequency is less straightforward than often assumed in the literature. The effect of 
an outer disc edge is to flatten the spectrum below the magnetic alignment frequency 
there. We also find a correlation between the variability amplitude and luminosity, 
similar to that found in some AGN. 

Key words: black hole physics - X-rays: binaries, galaxies: jets 



of the interest in scaling between stellar and supermassive 
objects. 

We stress, however, that the processes we discuss in this 
paper are fully relevant not only to binary X-ray sources and 
AGN, but to all objects containing accretion discs, such as 
cataclysmic variables and protostellar objects. 

Two main diagnostics are relevant for the variability of 
these sources. These are the shape of the power spectrum 
and the amplitude of the variability. Of course, the latter 
can be deduced from the normalisation of the former (e.g. 
Miyamoto et al. 1992) and so strictly speaking, these are 
just two aspects of the same diagnostic. 

The best observed black-hole candidate X-ray binary 
(BHC) is Cygnus X-l. In the low/hard state the normalised 
rms variability, a, is of order a few tens of per cent and the 
power spectrum can be approximated as a broken power 
law of the form p(f) oc f~ a where a ~ for 10 -3 Hz 
< / < 0.1 Hz, a ~ 1 for 0.1 Hz < / < f bl ~ 3 Hz, and 
a ~ 2 for / br < / < 10 3 Hz (Nowak et al. 1999; Revnivtsev, 



Gilfanov & Churazov 2000; Churazov, Gilfanov & Revnivt- 
sev 2001). This implies that the flickering variability occurs 
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predominantly on timescales in the range 0.3 s - 10 s. In 
the high/soft state the normalised rms variability, a, is of 
the order of a few per cent, and the spectrum is a similar, 
broken power law, except that now /b r — 10 Hz, and a ~ 1 
for / < /br- In this state, therefore, the flickering occurs pre- 
dominantly on timescales > 10 s. Other BHCs show similar 
behaviour. The power spectra of GX339-4 and GS2033+338 
in the low/hard state are given by Miyamoto et al. (1992) 
who note their strong similarity to Cyg X-l. The power 
spectrum of LMC X-l is shown by Nowak et al. (2001) to 
be similar to Cyg X-l in the high/soft state, and to have a 
normalised rms variability of a ~ 7 per cent. 

X-ray power spectra of Seyfert 1 galaxies also typically 
display a broken power law, with the break occurring at 
frequency /b r — 10~ 6 Hz, corresponding to a timescale of the 
order of a few days (Markowitz et al. 2003). For / < /b r the 
slope of the spectrum has typically a ~ 1 and for / > /br the 
slope of the spectrum has typically a ~ 2 (see also Vaughan 
& Fabian 2003; Vaughan et al. 2003). The normalised rms 
variability amplitudes range typically up to a ~ 10 per cent 
(Markowitz et al. 2003). 

Uttley and McHardy (2001) have drawn attention to 
a linear correlation between the actual rms variability, a, 
and the flux, F, in both BHCs and Seyfert 1 galaxies, im- 
plying that the normalised rms variability a — a/F is ap- 
proximately independent of source brightness, at least over 
a small range in flux. From the data for the BHCs Cyg X- 
1 and SAX J1808.4-3658, given in the paper by Uttley & 
McHardy (2001), we find that the relationship is approxi- 
mately of the form a oc F q where q ~ 0.3. 

The earliest models for the variability were in terms of 
shot noise (Terrell 1972); the light curves were reconstructed 
from a series of independent overlapping shots. Such a pic- 
ture, or variations thereof (e.g. Lehto 1989), can be made 
to model the data reasonably well, in particular the shape 
of the power spectrum. However, this scenario still lacks 
a physical picture for the variability. Since then, various 
ideas have been put forward, including accretion disc tur- 
bulence (Nowak & Wagoner, 1995), magnetic flares in the 
disc corona (e.g. Poutanen & Fabian, 1999), and MHD in- 
stabilities in the plunging region of the inner disc (Hawley 
& Krolik, 2001). These models all have a common prob- 
lem in accounting simultaneously for the range of timescales 
and the normalised rms variability amplitudes. The latter 
are typically around 10 per cent, and imply emission from 
the inner few disc annuli, where most of the accretion en- 
ergy is released (see for example the discussion by Bruch, 
1992; Fritz & Bruch, 1998). But these models also imply 
that characteristic variability takes place on local dynami- 
cal timescales. At the inner radii of accretion discs around 
~ 10A/q black holes, these relevant timescales are far too 
short, typically of the order of milliseconds. 

An insightful paper by Lyubarskii (1997) offers a way 
out of these difficulties. Lyubarskii considers small ampli- 
tude local fluctuations in the accretion rate at each radius, 
caused by small amplitude variations in the viscosity, and 
then considers the effect of these fluctuations on the accre- 
tion rate at the inner disc edge. A linear calculation shows 
that if the characteristic timescale of the viscosity variations 
is everywhere comparable to the viscous (inflow) timescale, 
and if the amplitude of the variations is independent of ra- 
dius, then the power spectrum of luminosity fluctuations has 



the form p(f) oc 1//. If the amplitude of the variations in- 
creases with radius, the slope of the power spectrum of the 
luminosity variations is steeper than —1. Lyubarskii (1997) 
notes that he has no physical model for the cause of such 
fluctuations. In particular, although the obvious candidate 
cause is the magnetic dynamo, the characteristic timescales 
for the dynamo are much shorter than the local viscous 
timescale. 

In this paper we use the idea (Livio, Pringle & King 
2003) that although the local dynamo timescale at each 
point in the disc is indeed short, the timescale for the dy- 
namo processes in a sufficient number of independent neigh- 
bouring disc annuli to produce a poloidal magnetic field co- 
herent enough to affect the accretion rate, by generating a 
disc wind or jet, can be quite long. That a disc wind/jet 
can be a significant driver of accretion in such objects is in 
line with the recent ideas about 'jet-dominated' states in 
black hole candidate binaries proposed by Fender, Gallo & 
Jonker (2003). We propose a specific, but highly simplified, 
model of an accretion disc. A magnetic dynamo process gen- 
erates the viscosity, but also, from time to time, produces 
a sufficiently well-ordered poloidal field. Therefore, the an- 
gular momentum is transferred outwards by two separate 
processes — the usual magnetic viscosity and a disc wind 
or jet. Numerical simulation of all the physical processes in- 
volved in such an intricate system is, of course, unrealistic 
at the current time, and we make no attempt to describe 
the full MHD problem in its three-dimensional complexity 
(e.g. Hujeirat, Camenzind & Livio 2002). Instead, we make 
a number of simplifying assumptions about the physical pro- 
cesses involved, and use these to construct a numerical de- 
scription of the time-dependence of such an accretion disc. 
We show below that this straightforward approach, based on 
simple physical ideas, can describe the flickering behaviour 
observed in accreting black hole systems. 



2 THE EQUATIONS FOR DISC EVOLUTION 

In this Section we derive and explain the equations we use to 
describe the time evolution of the disc. We simplify the prob- 
lem by adopting a one-dimensional formulation and using 
current ideas and simulations to introduce some of the basic 
physics. We measure disc radii in units of the inner radius 
-Rin, and assume an outer radius i? out large enough (typi- 
cally 20, 000-Rin) that the outer disc can act as a reservoir, 
allowing the inner disc to approach a quasi-steady state. 
We avoid undue complexity in these preliminary computa- 
tions by ignoring the details of internal structure and ther- 
mal effects within the disc. Instead, we assume a fixed disc 
opening angle, with thickness H oc disc radius R. We typi- 
cally take H/R — 0.08 for numerical convenience: a smaller 
and perhaps more realistic value would require longer run 
times, without changing the character of our results. We use 
a fixed grid with logarithmic spacing. We choose the zone 
width DR ~ H and thus TV ~ 130 zones. This is not strictly 
necessary, but gives sufficient accuracy for our current pur- 
poses — we note that the main errors in our treatment are 
likely to come more from our assumptions than from nu- 
merical inaccuracy. Indeed, using a formal numerical radial 
resolution which is finer than the disc thickness, is not likely 
to have physical meaning. In addition, the convenience of 
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the choice DR ~ H becomes apparent in Section \l. 2. II The 
disc material is assumed to have Keplerian angular velocity 
fi = {GM/R 3 ) 1/2 about a central point mass M, and we use 
units with GM = 1. In the present paper we consider a disc 
extending inwards close to the last stable orbit around the 
black hole, as is generally thought to occur in the high/soft 
X-ray state. 

2.1 The surface density 

We assume that the surface density E evolves because of 
viscous angular momentum transfer within the disc and be- 
cause of angular momentum loss in the magnetic wind or 
jet. Note that we use the terms jet and wind simply to mean 
a mechanism by which the disc can lose energy and angu- 
lar momentum, without being more specific. Thus we write 
(Pringle, 1981; Livio & Pringle 1992) 



as 
dt 



3.JL 
RdR 



(i) 



where Ur < is the radial velocity induced by angular 
momentum loss in the jet/wind (see Section 12.2.21 . Note 
that we assume that although the jet /wind removes angular 
momentum from the disc, it does not remove a significant 
amount of material. It is clear from global energy consider- 
ations that the wind cannot remove all of the disc material 
(put simply, all of the accretion energy cannot go into driv- 
ing all of the accreting material back to infinity). Hence, at 
the level of accuracy to which we are working, we neglect 
mass-loss to the wind in the surface density evolution. 

The viscosity v is assumed to result from magnetic 
torques generated by dynamo activity within the disc. We 
use here the standard Shakura-Sunyaev (1973) parametrisa- 
tion of that viscosity using the dimensionless parameter a 
in the form 



ac a H, 



(2) 



where c s is the (appropriately averaged) sound speed in the 
disc. Using hydrostatic balance we have approximately 



c s = HQ, 
and hence 
v = aH 2 n. 



(3) 



(4) 



Thus if there is no loss of angular momentum to the wind, 
(i.e. Ur = 0), then the evolution equation is a linear diffu- 
sion equation for the surface density E. We also note that 
by using this simple prescription and ignoring the thermal 
structure of the disc, we have for the time being put aside 
discussion of the role that hysteresis loops in the (y, E)-plane 
might play. These are the well-known S-curves which are 
thought to give rise to the outbursts of cataclysmic vari- 
ables and X-ray transients, and might well play a role in 
such systems (e.g. Belloni et al. 1997; Watarai & Mineshige, 
2003). In this paper we consider disc instabilities caused by 
the interaction of viscous and wind torques alone. 

We take a zero torque condition E = at the inner 
boundary R = 1, allowing accretion to occur freely through 
the inner boundary. At the outer boundary we impose a 
zero radial velocity condition. We use a simple first order 
explicit numerical scheme, using upwind derivatives for the 
advective Ur term. The scheme is written so that, except 



at the inner boundary, disc mass is conserved to machine 
accuracy. 



2.2 The magnetic field 

At each point in the disc we postulate a poloidal (vertical) 
field component B z generated by the disc dynamo produc- 
ing the viscosity v. We assume that this field drives the 
magnetic wind torques causing the inflow velocity Ur. We 
assume that the field is advected by Ur, and that it is also 
able to diffuse through the disc because of the magnetic dif- 
fusivity rj* resulting from the dynamo process. We discuss 
each of these assumptions below when describing their im- 
plementation in our disc evolution scheme. They imply an 
evolution equation for B z in the form 



OB 1 1 

-dT = ~RdR iRB * UR) + RdR 



dB z 
OR 



+ S B - 



(•») 



We should note here that our scheme neglects the advec- 
tion of magnetic field by the velocity induced by the viscos- 
ity v. The dynamo process generating B z in random fashion 
(Section |2~2~TI is an integral part of the process generat- 
ing viscous torques, and hence viscously driven accretion. 
This allows diffusion of B z radially through the disc (Sec- 
tion ^. 2.3"! . But, for unit Prandtl number, the timescales on 
which B z varies locally due to the dynamo (fc,jf2 -1 , where 
physically we expect kd ~ 10, Section r2.2.1|l and due to diffu- 
sion between neighbouring disc annuli (~ H 2 /rj ~ a~ 1 Q~ 1 , 
Section ^. 2. 3H are typically much shorter than the radial vis- 
cous diffusion timescale (r v i S c ~ R 2 jv). 

At the outer disc boundary we take B z = 0, allowing 
magnetic flux to diffuse outwards, but preventing inward 
flux advection. At the inner disc radius we assume zero flux 
loss through either diffusion or advection. The scheme is 
written so that poloidal flux is conserved to machine accu- 
racy through the disc, accounting for diffusion through the 
outer boundary and the dynamo source term. 



2.2.1 The source S B of B z 

We assume that the source of the local vertical field com- 
ponent B z is the disc dynamo process driving the viscos- 
ity. Specifically we assume that, in this respect, each annu- 
lus of the disc of width DR » H acts independently of all 
other annuli. We further assume that each annulus gener- 
ates large enough fields £?disc internal to the disc to produce 
the pseudo-viscous magnetic torques contained in the a— 
parameter. That is, we assume (Shakura & Sunyaev 1973) 



J>2 

-"disc 



(E/2H) c 2 



(6) 



and also that the dynamo process gives rise to a local vertical 
field B z of magnitude <£?disc. We assume, that for each 
annulus, B z can be modeled in terms of a stationary series 
resulting from a stochastic process. We model this process in 
terms of the Markoff scheme (e.g. Kendall 1976). This is the 
simplest linear autoregressive scheme other than a purely 
random series. In this scheme the n-th member of the series 
u n is given in terms of the previous by the relation 



-Qiu„_i + e„, 



(7) 
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where ot\ is a parameter, which for stability has modulus 
less than unity, and e n is a random variable with zero mean. 

This scheme produces a Markoff series with oscillations 
of a more or less regular kind. The mean number of series 
steps between oscillation peaks is 27r/cos~ 1 [— (1 + ai)/2], 
and the amplitude of the oscillations is such that 

varw = (°) 



In our simulations we choose ot\ = —0.5, so we expect 
3.45 series steps between peaks. We envisage the dynamo 
process to have some local canonical timescale Td(R). We 
assume that Td(R) = kdfl^ 1 with kd ~ 10 (e.g. Tout & 
Pringle 1992; Stone et al. 1996). We therefore apply the 
Markoff process through the source term S with the appro- 
priate frequency at each radius so as to give rise to that 
mean timescale. 

We vary the amplitude of the process by adjusting the 
variance of the term e n . To do this we use a dimensionless 
parameter B s . Then e is chosen to be a random variable in 
the range [B max , -B max ]. The value of -B max is proportional 
to B s and chosen so that if B s = 1 the resulting B z would, 
on average, produce a magnetic torque (see 12. 2. 211 driving 
an inflow velocity similar to that generated by the viscosity, 
i.e. v /R. Thus we take 



Bm»T = -:B a 7raE — 



H GM 



R R 2 



2.2.2 The jet torque and Ur 



(9) 



We assume that the poloidal magnetic field generates a 
jet/wind removing angular momentum and energy, but neg- 
ligible mass, from the disc, and so produces inflow velocity 
Ur. One way to proceed (cf. Lovelace, Romanova & New- 
man 1994; Lovelace, Newman & Romanova 1997) would be 
to assume that the inflow velocity in the disc is simply pro- 
portional to B 2 . However, since the source of the poloidal 
field is the dynamo process which has a radial scale in the 
disc of order ~ H <C R we feel that this is not physically 
justifiable. Instead, we follow the ideas of Tout & Pringle 
(1996; see also Livio 1996). They suggest that an inverse 
cascade driven by disc dynamics and reconnection can trans- 
form these small-scale poloidal fields into a global poloidal 
field smaller in magnitude but larger in scale. The compo- 
nent of the poloidal field with a scale of order ~ R could 
then produce a magnetically dominated outflow. 

To model these ideas in the current simulation we define 
two quantities 



and 



Sn + -Z> z R d R 



f 

J R 



RA 



R+A/2 B 2 RdR 



-A/2 



RA 



(10) 



(11) 



where typically we take the range A ~ R. Then the quantity 



Q(R) = (B z f/(Bl 



(12) 



is a dimensionless measure of the local coherence of the 
dynamo-generated field, and is such that < Q < 1. 

We could now simply allow the local magnetic torque 
to depend on (B z ) rather than B z (cf. also Livio & Pringle, 



1992) and assume that it gives rise to an inflow velocity, say 
Vr, given by 

(B Z ) 2 R 3 / 2 



Vr = 



(13) 



(GM) 1 /2 7r E' 

However, we must also consider the ability of a poloidal 
field to launch a wind or jet. Although a sufficiently large- 
scale poloidal field is probably able to give rise to a wind 
of some sort, Blandford & Payne (1982) point out that 
the wind is considerably enhanced if the poloidal field lines 
make a large angle with the vertical, i.e. if Br/B z is large 
enough (see also Ogilvie & Livio 2001). Such a field con- 
figuration results if the inflow velocity is comparable to the 
speed, ~ rj* / R, with which vertical field can diffuse outwards 
through the disc. Here rf is the effective magnetic diffusiv- 
ity in the disc, defined below in Section f2.2.3l This effect is 
discussed by Lubow, Papaloizou & Pringle (1994b) and by 
Lovelace, Romanova & Newman (1994) who conclude that 
it can drive strongly time-dependent disc behaviour. These 
ideas are confirmed by the 'accretion-ejection' events seen 
in the numerical calculations of Casse & Keppens (2002). 
Lubow et al. (1994b) suggest that the strength of the wind 
is a sensitive function of the quantity X = RVr/t/* . Lubow 
et al. (1994b) argue that for a locally isothermal disc the 
loss from the disc is oc exp(— 1/X 2 ), and that for large X 
the loss timescale is the local dynamical one. In the simula- 
tions reported here we do not use the non-analytic function 
exp(— 1/X 2 ) (sometimes called the April Fool Function, be- 
cause, although all its derivatives exist at the origin, it is not 
equal to its null Taylor series about that point). Instead, for 
numerical reasons we approximate the sensitivity to X by 
assuming that the radial velocity driven by the jet/wind 
obeys 



Ur = Vr{1 + {/lpp - 1} 



X 2 



X 2 + 4 



(14) 



Here the parameter /lpp, assumed > 1, represents the 
maximum enhancement to the torque that this effect can 
achieve. This formula has the behaviour that as X — > oo, 
Ur —> ,/lppVr, and Ur tends to Vr at small X. Note that 
setting /lpp = 1 removes this effect. 



2.2.3 The magnetic diffusivity rj* 

We assume that the disc dynamo producing the effective 
viscosity v also gives rise to an effective magnetic diffusivity 
rj. It seems reasonable to assume that the Prandtl number is 
of order unity, and thus that rj — v. However, if the poloidal 
field is aligned over a scale of order the radius R then it 
diffuses through the disc at a rate enhanced roughly by a 
factor of R/H (van Ballegooijen 1989; Lubow, Papaloizou 
& Pringle 1994a). Of course, the actual rate at which the 
poloidal field diffuses through the disc depends on the details 
of the configuration of the field which we do not attempt to 
model here. We have argued above that the quantity Q is 
a measure of the local coherence of the dynamo-generated 
field. Thus we approximate these effects by assuming that 
B z is subject to an effective magnetic diffusivity rf acting 
in the radial direction, where 



rj* = 77max(l, QR/H), 

and Q is defined in equation 11211 . 



(15) 
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3 RESULTS 

Even though we have approximated much of the essential 
physics of a full MHD disc calculation for reasons of simplic- 
ity, we nevertheless require a number of parameters to spec- 
ify our results. These are a, H/R, B s , /lpp, k d — T d (R)£l(R), 
and Rout/Rin- The resulting behaviour of what is a highly 
nonlinear system is sufficiently complicated that we do not 
attempt, in this paper, to undertake a full survey of param- 
eter space. Rather we content ourselves with fixing most of 
the parameters. We take Rout/Rin = 20,000 and note that 
the actual value of this is unlikely to play a significant role 
provided that it is 3> 1. For the runs we discuss below (un- 
less indicated otherwise) we take B a = 10, since this value 
seems to ensure that typical values of the mean poloidal 
field, {B z ), are large enough to have a noticeable effect on 
the angular momentum loss, but small enough (factors of a 
few times less than the disc field Bdisc) to be physically plau- 
sible. Similarly, unless stated otherwise, we take /lpp = 3, 
since this value ensures that the physical effect of the jet 
torque plays some role without being unphysically domi- 
nant. We also choose parameters such that H/R ~ 0.08 (see 
below). We then experiment with varying the remaining pa- 
rameters a and kd — Td(R)£l(R). 

In each run, we start from an initial condition in which 
the disc surface density is constant at all radii and equal 
to the arbitrary value of unity (E = 1). [Note that we can 
do this because v is independent of surface density, hence 
the evolution equation Q for E is linear in the absence of 
wind/jet angular momentum.] We start with vertical mag- 
netic field B z = everywhere, but then take one initial step 
in the Markoff process at each annulus. 

In order to get some idea of observable disc parameters 
we compute estimates of what we denote as a 'disc lumi- 
nosity' Ld and a 'wind/jet luminosity' Ljet, as functions of 
time. We regard the disc luminosity as some measure of the 
radiative emission from the accretion flow. However we do 
not model the conversion of gravitational energy into radi- 
ation, and the consequent division of this energy between 
hard and soft X-ray components. It may be that the instan- 
taneous wind/jet luminosity is generally unobserved, except 
perhaps when rapid variations give rise to internal shocks in 
the outflow (see, however, Malzac 2003). 

We take the disc luminosity as caused solely by viscous 
dissipation in the disc, and hence (e.g. Pringle, 1981) take 
each disc annulus DR to contribute an amount 

SL d = -v^^DR. (16) 
4 R? 

We estimate the rate at which energy is removed from 
each annulus by the wind/jet from the disc by assuming 

CM 

8L jct = (—Ur)tt'E——DR. (17) 
H 

Note that these units are arbitrary in the sense that 
the units for E are arbitrary, but that the relative value of 
the two quantities does have meaning. The total disc and 
jet/wind luminosities are then obtained by summing these 
over all the disc annuli, although evidently most of the con- 
tribution comes from the innermost radii. In the figures be- 
low we present values for these quantities averaged over 10 4 
code time units (~ 10 s for a central black hole of 14 Mq). 
We note that for our canonical parameter values of B a — 10 



and /lpp = 3, the disc and jet luminosities are approxi- 
mately scaled versions of each other. 

We emphasize that we have implicitly ignored the mag- 
netic energy represented by B z , which is created/destroyed 
by the dynamo process, and is dissipated by the diffusivity 
rj. One reason for doing this is that we have no proper model 
for the processes that generate and destroy B z or that de- 
termine its physical scale. However, even when Bdisc ~ B z , 
the local energy flow (per unit disc area) into and out of the 
poloidal field component is of order Q mag ~ H (B 2 / 4tt) / r d 
and is therefore less than the rate at which energy is dis- 
sipated in the disc by viscous effects (which depends on 
a oc B% iBC ) by a factor of r^fi = kd S> 1. Thus, given the 
simplifications we have made, neglect of the poloidal energy 
is reasonable at this level to a first approximation. 

We take the unit timescale for the computations so that 
the orbital timescale at the inner edge (R — 1) is 2n. In or- 
der to make comparisons with observations more obvious, 
in the diagrams and in the discussion below we have taken 
parameters such that M = 14M©, and Ri n = 6GM/c 2 
so that Ri n = 124(M/14M ) km and the time unit is 
1.016(A//14Mq) ms. Runs are carried out typically for a 
time of 1.6 x 10 7 (M/14M Q ) ms. 

3.1 Timescales and frequencies 

Before presenting the results of some of our simulations 
we briefly discuss some relevant timescales. Thus, we con- 
vert code time units to milliseconds. For central masses 
other than 14 Mq the inverse angular frequency varies as 
Q- 1 oc M~ 1/2 R 3/2 , thus from 1O~ 3 (M/14M ) s at the in- 
side, R = 1, to approximately 2.8 x 10 3 (M/14M Q ) s at 
R = 20,000. Similarly, the dynamo timescale Td, taken 
equal to fc^fi -1 , varies as Td oc Af _1//2 J? 3/2 from a value 
of 0.01(fc d /10)(M/14M Q ) s at R = 1 to a value of 2.8 x 
10 4 (fc d /10)(M/14M Q ) s at R = 20,000. We define the 
dynamo frequency as fd = t// 1 . The maximum value of 
this frequency occurs at the inner edge and is equal to 
/^ ax = 100(fc (i /10)- 1 (M/14M Q )- 1 Hz. 

The viscous timescale t v is given (e.g. Pringle 1981) by 

t v « ^{R/Hfar 1 . (18) 

Thus for the value of H/R « 0.08 (independent of radius) 
adopted in these calculations, the viscous timescale varies 
as t„ oc M~ 1/2 B? f2 oT l from 16(M/14M e )/(a/0.01) s at 
R = 1, to 4.5 x lO 7 (Af/14M )/(a/O.Ol) s at R = 20,000. 
For the calculations presented here this means that at times 
t > 1, 000 s, the central regions of the disc (say, R < 20) where 
almost all the energy is generated, and where all the "action" 
is taking place, act as if they were being fed at a steady rate 
by a disc with unit surface density. The timescale on which 
magnetic diffusion works is at most the same as the viscous 
one, and can be less by a factor of up to R/H 12.6. With 
the viscous timescale we can associate a viscous frequency 
f v = t~ . The highest value of the viscous frequency occurs 
at the inner edge and for the disc parameters we use here is 
given by /™ ax = O.O625(a/O.Ol)(A//14M )- 1 Hz. 

The final relevant timescale T mag is that for a sufficient 
number of neighbouring annuli to generate an aligned mag- 
netic field, that significant local energy loss to the jet/wind 
occurs. This timescale is discussed in Livio, Pringle & King 
(2003). From our assumptions, the number of zones needed 
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Figure 1. The left panel (a) shows representative jet (top) and disc (bottom) lightcurves for kd = 10, a = 0.06. The lightcurves are 
normalised by the mean luminosity in the segment shown, and binned to 10s time resolution. The right-hand panel (b) shows the disc 
power spectrum calculated from a 16000s lightcurve. In this and all subsequent figures we give la error bars. < Lj et > / < idisc >= 1-37 
for this case. 
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Figure 2. Lightcurves (a, left) and power spectrum (b, right) for kd = 10, a = 0.006. Lightcurves have been normalised and binned as 
per the description of Figure 1. < Lj ct > / < ^disc >= 3.38 for this case. 



at radius R corresponds to the number of neighbouring an- 
nuli, width ~ H, required to span a radial scale of order R. 
Let this number be N + 1. Then, assuming for simplicity 
that B z can take only the values ±1, if each zone acts inde- 
pendently, the probability of all TV + 1 zones being aligned 
is 1/2 N . To generate a timescale for this we assume that all 
the zones at radius R are randomly realigned every Td(R)- 

To give specific examples, for most of the simulations 



we take the number of grid zones to be N sf id = 130 which 
gives H/R = 0.079 ^ 0.08 and N = 8. Thus the timescale 
on which the innermost 9 zones become aligned is approxi- 



mately 



2.56(fcd/10) s. Since r mag cx i? 3//2 the process 



occurs more slowly at larger radii. Thus we expect the pro- 
cess to be dominated by timescales at the smallest radii, 
but modulated by longer timescales corresponding to larger 
radii. The associated magnetic alignment frequency is given 
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by /„ 



As for the other frequencies the highest 



magnetic frequency occurs at the inner edge and is given by 
C = O-mkd/wy^M/UMe) Hz. 

As noted by Livio, Pringle & King (2003), the quan- 
tity T m ag depends sensitively on H/R. For example, if we 
consider H/R = 0.135 we find N — 4, and r mag is reduced 
by a factor of 16. However in the numerical method of this 
paper we use the fact that DR — H. Thus, changing H/R 
would require a different numerical grid, and also change 
the viscous timescale (for the same value of a). To investi- 
gate the effects of changing the magnetic timescale, without 
altering other parameters, we here simply vary the parame- 
ter kd- Of course in reality kd is determined by the physics 
of the disc. This procedure is designed to mimic the effect 
of changing H/R in a numerically convenient fashion. Ac- 
cordingly, we have carried out runs spanning the parameter 
space 0.001 < a < 0.6 and 10 < k d < 1000. 

3.2 Power spectra 

As we noted above, if kd = i~dQ is independent of radius, and 
if, as we have assumed, H/R and a are also independent of 
radius, then the ratio of the two relevant timescales r mag and 
t v is also independent of radius. Thus we expect our results, 
appropriately scaled, just to depend on the dimensionless 
ratio 



2 R ' H k d a 
{R/HY 



(19) 

T d 2 R,H (Livio, Pringle 



where we have used 11811 and r mag 
& King 2003). 

The power spectrum calculated for each lightcurve is 
normalised using the prescription of Miyamoto et al. (1991), 
namely 



215, 



n 2 T 



(20) 



where Sj is the Fourier amplitude at frequency fj, 1Z is the 
mean luminosity of the lightcurve, and T the total duration 
in seconds. We have applied a logarithmic frequency binning 
to all the PSDs, with A/// = 0.1. 

3. 2. 1 7" m ag — Tv 

In Figure 1(a) we show the light curves for Ldisc and Lj e t and 
in Figure 1(b) we show the normalised power spectrum of the 
disc luminosity for a = 0.06 and k d = 10. For these values 
of a and k d , r mag (i? = 1) ~ r„(R = 1) ~ 2.7/(M/14M ) s, 
or, equivalents, ^ /T* ^ 0.4(M/14M Q )" 1 Hz. From 

Figure 1(b) we see that at frequencies / < /™ ag ~ /™ ax 
the power spectrum has the form p(f) oc 1//. This corre- 
sponds exactly to the situation described and predicted by 
Lyubarskii (1997). At each radius the local inflow rate is per- 
turbed randomly with a characteristic timescale which is ap- 
proximately equal to the local viscous timescale. These per- 
turbations then flow inwards on the local viscous timescale, 
and when they reach the inner few radii contribute to the 
total disc luminosity with a fluctuation on that timescale. 
The fact that the 1/f spectrum extends to frequencies as 
low as 10 _3 /™ ax implies that contributions to the variabil- 
ity come from radii as large as R ~ 100. In Lyubarskii's 
(1997) model, the power spectrum drops exponentially at 



frequencies above /™ ax because his assumed perturbations 
have no power at frequencies higher than this. He reached 
this conclusion because, as he remarked, he had no specific 
physical model for the perturbations. In the model we con- 
sider here there are perturbations to the local inflow rate 
at all timescales down to the local dynamo timescale at the 
inner disc edge, Td(R ~ 1), corresponding to a frequency of 
fd — 100 Hz. At frequencies / > /™ ax the only perturbations 
affecting the disc luminosity are those which take place at 
radii close to the inner edge where the potential well is deep- 
est. Thus the shape of the power spectrum in this frequency 
regime depends simply on the power spectrum of the applied 
random perturbations. From Figure 1(b) it is apparent that 
in the model presented here the power spectrum steepens for 
frequencies above about /br ~ 0.2 Hz sa /m ag . The spectrum 
has a slope of d logp/d log / = -2 at / « 2/(M/14M ) Hz, 
and = -3 at / « 2O/(M/14M ) Hz. We see a similar be- 
haviour, suitably scaled in frequency space, for the run with 
a = 0.006 and k d = 100. 



3.2.2 



< T„ 



We show in Figures 2(a) and 2(b) the light curves and power 
spectra for a similar model with k d = 10, but now with a = 
0.006. Thus in this case we have f£g = O.4/(M/14M ) Hz 
> /™ ax = 0.04/(M/14M Q ) Hz. As mentioned above, we 
have the ratio /""//msg = 0.1 independent of radius. In 
this case the power spectrum is nowhere of the form p(f) oc 
1/f, and indeed is quite well represented by p(f) cx l// 2 
at all frequencies. We see similar behaviour for the case 
k d = 10, a = 0.001, for which = 0.017. In this 

case the magnetic torque fluctuations throughout the disc 
occur locally much faster than the timescale on which vis- 
cous forces can react and inflow can occur. This implies that 
the magnetically induced variability has only a local effect. 
Thus the disc luminosity, which is dominated by contribu- 
tions from the inner regions, is only affected by magnetic 
torque fluctuations in those inner regions. 



3.2.3 



> T„ 



In contrast with the above, in Figures 3(a) and 3(b) we 
show the light curves and power spectra for a similar model, 
with kd ~ 10, except that a = 0.6. Thus in this case we 
still have /™ a g = 0.4/(M/14M Q ) Hz, but now /™ ax = 
4/(M/14M Q ) Hz > /™ ag - In this regime the power spec- 
trum is approximately of the form p(f) oc 1/f for fre- 
quencies /</br, where /b r ~ 2/™ ag , independent of the 
value of /™ ax . In the frequency range f br < f < /™ ax , the 
power spectrum is roughly of the form p(f) oc l// 2 , and the 
power spectrum falls off more steeply at higher frequencies 
/ > ,f™ ax - These findings are also illustrated by the cases 
a = 0.06, k d = 100 and a = 0.06, k d = 1000 (Figure 4). In 
this case, because the viscous timescale is short, all magnetic 
torque fluctuations can be communicated to small radii, and 
so can contribute efficiently to fluctuations in the disc lumi- 
nosity. At frequencies lower than about /°™ the situation 
is still essentially as described by Lyubarskii (1997) and the 
power spectrum is of the form p(f) oc 1/f. However, at 
frequencies higher than this, but still below Z™ 2 ", the fluc- 
tuations carried inwards to small radii have ever lower am- 
plitudes, so the power spectrum steepens. At frequencies 
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Figure 4. Power spectra for a = 0.06, implying /™ ax = 0.4 Hz, and the cases ka = 100 (left) with /, 
: 0.004 
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higher than Z™' 1 " however, there is no time for fluctuations 
to be carried inwards from larger radii, and all we see is the 
effects of local fluctuations occurring in the innermost disc 
radii. Hence the power spectrum steepens still further. 



3.2.4 Effect of the outer disc edge 

In the computations presented so far, all the processes have 
varied with radius in a self-similar manner. Thus, for exam- 
ple, the ratio /„//ma g was independent of radius. For this 
reason, the structure of the power spectra has been rela- 
tively straightforward to interpret. In general, of course, the 
situation is likely to be more complex, with, for example, 
both /„ and /mag being intricate, and different, functions of 
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Figure 6. The evolution of the disc and jet throughout a large amplitude flare (t<i = 100, a = 0.6). Four selected times are marked on 
the lightcurve. The top panel shows the disc lightcurve segment covering the flare event. The four lower panels show the disc surface 
brightness, surface density and jet luminosity for the four selected times. 
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Figure 5. The ratio < Lj ct > / < Ldisc > versus the viscosity 
parameter a, for kd = 10. 
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Figure 7. The power spectrum for the case a = 0.06, kd = 10, 
for which the magnetic wind activity has been turned off at radii 
R > 100, corresponding to a frequencies of / < f sa 0.0013 Hz. 



radius. In this Section we illustrate one simple complication, 
by setting the dynamo wind activity to zero outside some 
radius, _R m ax- This radius, can be viewed for example as 
the outer edge of the disc. Alternatively, it can be regarded 
as a radius outside which the magnetic alignment timescale, 
Tmag becomes very long, because the disc thickness suddenly 
drops outside that radius. 

In Figure|7|we show the power spectrum for a run which 
is identical to that shown in Figure 1 (with kd = 10 and 
a — 0.06), except that dynamo wind activity is turned off 
for radii R > J? ma x = 100. Since we have assumed that the 
presence of a wind/jet depends on magnetic alignment activ- 
ity over about a factor of two in radius, we should expect this 
to give rise to a feature in the power spectrum at a frequency 
related to f D « /„(-R max /2) = / mag (i? ma x/2) = 0.0013 Hz. 
As can be seen by comparing Figures 1 and t ne effect of 
truncating wind/jet activity at R = 100 is to cause a flat- 
tening of the power spectrum from P(f) tx 1/f to P(f) ~ 
const, at a frequency / ~ 0.002 Hz. 

3.3 Lightcurves 

In Figures 1 - 3, we have shown some typical samples of 
lightcurves for both disc and wind/jet emission for kd — 10 
and for a = 0.06, 0.006 and 0.6. From the Figures it is ev- 
ident that the disc and wind/jet luminosities follow each 
other fairly closely, except that the wind/jet emission is 
much less smooth. This comes about because the wind/jet 
emission depends directly on the local magnetic field config- 
uration, whereas the disc emission depends on the local disc 
mass flux which is also mediated by the viscosity. 

^From the simulations we find, as we increase a, that 
the amplitude of the flickering increases (see Figures 1-3). 
We discuss this in Section ^. 3. 2l below. We also find that the 
ratio of average jet luminosity to average disc luminosity de- 



creases (See Figure 5). As can be seen from Figure 5, for the 
parameters we have chosen (Lj et ) / (I/disc) is a slowly varying 
function of a and is unity for a m 0.1. Since the value of the 
poloidal field is scaled to ensure that the inflow caused by 
the wind/jet is on average similar to that caused by the vis- 
cosity, it is not immediately obvious what causes this effect. 
One possibility is that as a increases, keeping the dynamo 
timescale (measured by kd) the same, the wind /jet is able to 
sweep material, and magnetic flux, inwards more efficiently, 
and therefore to do so at lower values of the poloidal mag- 
netic field. 

3. 3. 1 Large flares 

In addition to the continual small-scale flickering shown by 
our simulated lightcurves, we also see occasional large ampli- 
tude flares. In Figure 6 we show an example of one of these, 
for a run with the parameters a — 0.6, kd = 100. Evidently 
this flare has a cause similar to that of the small amplitude 
flickering, except that it is instigated from a much larger 
radius (R ~ 100), is much more coherent, and so is of much 
larger amplitude. Such large flares, which start at large ra- 
dius, occur seldom because the timescales at large radii are 
correspondingly longer. The Figure shows the disc luminos- 
ity light curve (we note that the jet luminosity is similar but 
more spiky) , as well as the state of the disc at four particular 
epochs (marked in the luminosity plot). In these four plots, 
we note that the disc surface brightness and disc surface den- 
sity are simply proportional to each other, and are relatively 
smooth, whereas the jet surface brightness is much more er- 
ratic, reflecting the stochastic nature of the local poloidal 
magnetic field. Box 1 and Box 4 show the disc pre- and 
post-flare. Comparison of these two Boxes shows that the 
flare depletes the disc at radii R < 100. The timescale for the 
recovery of the disc luminosity post-flare indeed corresponds 
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Figure 8. Correlation of the rms variability with flux for three models: a = 0.006 (left), a = 0.06 (centre) and a = 0.6 (right), = 10 
in all cases. 



to the viscous timescale at a radius of R ~ 100. We note that 
observations of GRS 1915+105 (Belloni et al, 1997) show 
evidence for the viscous refilling of a disc depleted by flares. 
Box 2 shows the state of the disc at the start of the flare. 
The incipient dip in surface density at around R ~ 80 is the 
result of a coherent magnetic structure, seen in the jet sur- 
face brightness, starting an avalanche of material inwards. 
Box 3, close to the peak of the flare, shows that this wave 
of material has moved further inwards, and that the surface 
density in the inner disc has risen accordingly. We note that 
this behaviour demonstrates the ability of angular momen- 
tum losses in a disc wind/jet to give rise to a wave-like surge 
of inwardly propagating material, in contrast to the diffusive 
behaviour of the more usual transfer of angular momentum 
through viscosity. This can be seen clearly in Equation (1) 
which governs the evolution of surface density. There, the 
wind/jet term involving Ur contains a single radial deriva- 
tive and so is wave-like, or hyperbolic, in character, whereas 
the viscous term, involving v has a double radial deriva- 
tive and so is diffusive, or parabolic, in character. The idea 
that angular momentum loss to a wind might produce such 
avalanche/flaring behaviour has been discussed previously 
by Lubow, Papaloizou & Pringle (1994b), by Lovelace, Ro- 
manova & Newman (1994), and by Livio, Pringle & King 
(2003). 



3.3.2 Dependence of flickering amplitude on luminosity 

Uttley & McHardy (2001) have shown that in some AGN 
there is a linear correlation between the amplitude of the 
flickering and the X-ray flux. To test whether similar be- 
haviour is evident in our model, we divide our lightcurves 
into 100 second segments and calculate the mean flux (Fj) 
and rms flux (<7j) for the jih segment as follows: 

P i = JfT, F *<> ( 21 ) 

i 

and 



where the index i runs over all time bins in the segment 
j. The time bins in our lightcurves are 0.1s duration, and 
thus the o~j we calculate can be roughly interpreted as the 
rms variability in the frequency range 0.01 - 10 Hz. The 
normalised rms variability is then defined as a — a/F. 

In Figure 8 we show plots of a as a function of disc lu- 
minosity for the three models whose lightcurves we show in 
Figures 1-3. These have kd = 10 and have a = 0.006, 0.06 
and 0.6. In each plot the variability amplitude correlates 
approximately in a linear fashion with the luminosity. The 
normalised rms variability is then measured by the gradient 
of this relationship and is therefore roughly a constant for 
each value of a. Thus, we find that the normalised rms vari- 
ablity is 3.5 per cent for a = 0.006, 14 per cent for a = 0.06, 
and 28 per cent for a = 0.6. This increase of normalised rms 
variability with a is also evident from the power spectra 
shown in Figures 1 - 3. As was apparent in our discussion 
of these power spectra, an increase in a leads to an increase 
in the efficiency with which fluctuations at large radius in 
the disc are swept inwards to small radii, and thus to an 
increase in the amplitude of the variability of the disc lumi- 
nosity, which is produced mainly at small radii. 



4 DISCUSSION 

Although the disc behaviour we investigate here is relatively 
complex, one needs to remember that even within the limits 
of the simple model we have introduced, there are a num- 
ber of physical processes which we have not addressed, but 
which are nevertheless likely to play an important role in 
any observed source. 

One major omission is any consideration of the thermal 
behaviour of the disc. This can have a major impact because 
the magnetic alignment timescale is a sensitive function of 
R/H, with r mag tx r d 2 R/H (Livio, Pringle & King 2003). 
This has a number of implications. First, since in most discs 
the ratio H/R is not a constant, outflow activity is likely 
to be strong function of radius — for example close to the 
Eddington limit, the inner parts of the Shakura-Sunyaev 
disc (1973) have H « const, and hence H/R oc 1/R, which 
would imply enhanced outflow activity in the central disc re- 
gions. Indeed, we already know that because jets tend to be 
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launched mainly from the centres of accretion discs, the out- 
flow behaviour is likely to be strongly non-self similar (See 
the discussion by Pringle 1993; Livio 1996; Price, Pringle & 
King 2003). Second, the ratio H/R is in general an increas- 
ing function of the accretion rate, or mean disc luminosity. 
Thus in any given source it is likely that outflow activity, 
and hence the nature and magnitude of the variability, will 
depend strongly on the brightness of the source. This effect 
comes on top of the one presented in Section [3.3.21 Indeed, 
the data presented by Uttley & McHardy (2001) show some 
indication that the relation is steeper than linear. Third, a 
local enhancement of inflow rate caused by outflow activ- 
ity also leads to enhanced local dissipation within the disc. 
This leads in turn to a local increase in disc thickness, a 
marked decrease in r mag and thus to enhanced outflow ac- 
tivity. Thus, it seems likely that local outflow activity is sig- 
nificantly self-enhancing. This phenomenon would lead to 
much more pronounced activity than we find in our current 
models. Fourth, the thermal hysterisis behaviour which can 
give rise to outbursting behaviour on its own (as proposed 
for example for GRS1915+105 by Belloni et al. 1997, and 
Watarai & Mineshige 2003) could be considerably affected 
by enhanced outflow behaviour particularly during the hot 
state when the disc thickness is increased (see also the dis- 
cussion by Livio, Pringle & King 2003). 

Apart from thermal considerations, there are also likely 
to be further complications which will only come to light 
when more detailed numerical simulations of the processes 
involved become possible. One possible effect here is the rela- 
tionship between the strength of the poloidal magnetic field 
and the dynamo process. The dynamo process is thought to 
be driven by the magneto-rotational instability (MRI; Bal- 
bus & Hawley, 1991). This instability shuts off if the poloidal 
field becomes too strong. Suppose then that the dynamo ac- 
tivity at large radius produces (as we have assumed here) 
a poloidal field sufficient to drive an outflow. This causes 
wave-like inflow in the disc, carrying with it the poloidal 
field. As this wave of material moves to smaller radii, the 
poloidal field becomes advected with it and compressed, this 
process being mitigated by the effective diffusivity. If at some 
radius the field becomes strong enough to shut off the MRI, 
then at that radius both the effective disc viscosity and the 
effective disc diffusivity are significantly reduced. At that 
point, the poloidal field is then trapped by the disc, and the 
wave-like inflow becomes an unstoppable runaway. If this 
kind of behaviour occurs, then the avalanche-like behaviour 
we have found in the above models could be a severe un- 
derestimate of the violence and frequency of major flaring 
behaviour. 



5 CONCLUSIONS 

In this paper we have proposed a simple model to describe 
how the local dynamo process in an accretion disc, which 
produces the effective disc viscosity, can also affect the disc 
evolution by driving angular momentum loss in the form of 
a wind or jet. While the disc viscosity produces the usual 
smooth inflow, the accretion rate driven by the wind/jet is 
highly stochastic. This comes about because of our assump- 
tion that efficient angular momentum loss to an outflow can 
only occur when the local poloidal field produced by the disc 



dynamo is sufficiently coherent. We assume that each disc 
annulus (of size comparable to the local disc thickness H) 
produces a poloidal field whose size and direction varies as 
a random walk with some characteristic timescale (Markoff 
process). Only when a number ~ R/H of such neighbouring 
annuli produce sufficiently aligned poloidal fields can out- 
flow and angular momentum loss occur. 

We have implemented these ideas in a simple set of 
equations that determine the evolution of the disc, and 
hence the disc and jet(outflow) luminosities. We have chosen 
the parameters of the model (B a and /lpp) such that the 
poloidal field strength produced by the disc, while smaller 
than the tangled field within the body of the disc, is locally 
large enough to give a significant outflow. Thus, typically, 
in our models, Z/j e t and Ldi sc are comparable in magnitude. 
From an observational point of view this is not unreasonable 
(e.g. Fender et al. 2003; D'Elia, Padovani & Landt 2003). 
We have chosen simple parameters for our disc models, with 
constant ratio H/R and constant viscosity parameter a. In 
addition, we have assumed that local dynamo timescale is 
a constant x the local dynamical timescale f2 _1 in the disc. 
These assumptions imply that the disc is self-similar and 
that the local viscous timescale t v is proportional to the 
timescale r mag on which the local poloidal field alignment 
occurs. 

We find that in general the disc luminosity shows ir- 
regular behaviour. Small amplitude flickering is superposed 
on variations in the average luminosity, together with oc- 
casional large amplitude flares. We have investigated the 
properties of the flickering with the goal of obtaining an un- 
derstanding of what the flickering might be able to tell us 
about the properties of the disc. Lyubarskii (1997) showed 
that the power spectrum has the form p(f) oc 1// if the 
magnitude of viscosity varies on a timescale everywhere the 
same as the local viscous time (r v = l/f v )- In our case 
the first timescale is that for angular momentum loss in the 
outflow T mag = l//ma g . This is the case shown in Figure 
1, where /„ = /mag- The 1// spectrum continues up to a 
break frequency /b r corresponding to the inner disc edge, 
beyond which the spectrum steepens. For cases in which the 
local magnetic alignment timescale was longer than the lo- 
cal viscous time (i.e. / m a g < fu) we find a similar behaviour, 
with the frequency /t, r corresponding to the local magnetic 
alignment timescale at the inner disc edge. We do not agree 
with the assumption found in the literature (e.g. Vaughan 
& Fabian 2003) that the break frequency corresponds to the 
viscous timescale at a particular radius. If the local magnetic 
alignment timescale is everywhere shorter than the local vis- 
cous timescale, the spectrum is steeper than 1/f through- 
out (i.e. red noise). We have also shown that the effect of an 
outer disc edge is to produce a flattening of the spectrum 
for frequencies / < f , where f is the magnetic alignment 
frequency at that edge. 

We have investigated how the amplitude of the flick- 
ering varies with luminosity and with a. We find that in 
any one run, i.e. at fixed a, the normalised rms amplitude 
is approximately constant, in line with the findings of Utt- 
ley & McHardy (2001) from observation. We also show that 
the normalised amplitude of the flickering increases with in- 
creasing a, varying in our models from about 3 per cent for 
a = 0.006 to almost 30 per cent for a — 0.6. This proba- 
bly comes about because the larger the viscosity, the more 
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rapidly stochastic behaviour can be swept inwards from the 
outer disc regions. Thus, it may be that large amplitude 
flickering is a symptom of a large disc viscosity. 

We have also demonstrated the possibility of angular 
momentum loss to an outflow giving rise to avalanche-like 
behaviour resulting in a large flare in luminosity. This is be- 
cause angular momentum loss directly to an outflow changes 
the characteristic nature of the equation governing the evo- 
lution of the disc from diffusive to wave-like. We note that 
much of the flickering behaviour also comes about because 
of the wave-like nature of the response. This not only en- 
hances, but also keeps coherent, any perturbations to the 
disc. 

Finally, we stress again that, although we have applied 
our findings to accretion discs around black holes (BHCs 
and AGN), the whole range of phenomena found here, and 
much more (for example, interactions between the disc and 
the central object's surface and/or its magnetosphere) , are 
likely to occur in all systems in which accretion discs are 
found. 
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